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Abstract 


In this preliminary study of plasma effects on high-frequency electromagnetic 
wave propagation, the electromagnetic interaction with the shoc'k-induced plasma 
flow field generated around a l)lunt Innly capsule as it encounters the planetary 
atmosphere of Mars is examine<l. 

From two mathematical models of the Martian atmosphere (VM-4 and VM-8) 
and the entr>'- trajectory characteristics of the capsule, it appears that the communi- 
cation links to and from the capsule will be blacked out during entry because of 
the free electron concentration in the wake region of the capsuh-. 

Estimates of the free electron concentration in the wake region indicate peak 
values of approximately 8X10” e /cc for the \'M-4 model atmosphere and 
peak values of approximately 8 X 10'” e'/cc for the V^M-8 model atmosphere. 
Both of these values are above the critical electron exmeentrations of 6.53 X 10'“ 
e*/cc for S-bund transmission at 2.295 GHz and the critical electron concentration 
of 1.99 X 10“ e'/cc for transmission at 400 Mhz. Tliese frequencies are the trans- 
mission frequencies for communication links between the capsule and a relay bus 
and directly between the capsule and the earth, respectively. 
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A Martian Entry Propagation Study 


I. Introduction 

Part of the objectiv*. ' a mission to the planet Mars is 
to obtain data about th(* Nkartian atmosphere and surface 
conditions.’ To obtain these data, an instrumented capsule 
will be soft-landed on the surface of the planet. 

As the capsule enters the Martian atmosphere, a re- 
sistance* to its motion by the atmospheric gases will result 
in a transfer of energy from the capsule to tlie gases and 
a subsequent ionization of the gases. Tlie degree of ioni- 
zation will depend on the composition and properties of 
the Martian atmosphere and the entry-trajectory char- 
acteristics of the capsule. If the degree of ionization is 
substantial, a shock-induced envelope of ionized gases 
will form around the capsule as it desc*ends through the 
planetary atmosphere. 

Tlie envelope of ionized gases surrounding the capsule 
constitutes a plasma flow field, i.e., a macroscopically 
neutral ionized gas consisting principally of free electrons, 
free ions, and neutral particles. One effect of the plasma 
on electromagnetic wave propagation is an attenuation 


'“Mars *71 Technical Study.” Internal dcKument. 


of the wave, the .severity of the attenuation depending 
largely on the degree of ionization of the gases consti- 
tuting the plasma. At hypersonic entry velocities, the 
attenuation can be sufficient to cause total lo.vs of signal, 
or blachmt, even at extr<‘mely high frequencies. 

A preliminary examination of the various models of the 
Martian atmosphere (Ref. 1) and the entry-trajectory 
characteristics of the capsule* indicates that the signal 
attenuation will lie sufficient to cause blackout during 
some period of the entry pha.se of a Mars mission. The 
objective of this study is to determine the actual degree 
of signal attenuation or blackout to be expected during 
entry of a blunt capsule into the Martian atmosphere. 

A model to describe the actual plasma environment 
generated around the capsule as it enters the planetary 
atmosphere will be developed from two mathematical 
models of the Martian atmosphere (the VM-4 and V^M-8) 
that have been chosen to represent a reasonable range of 
expected conditions, and from the entry-trajectory char- 
acteristics of the capsule. With the use of the plasma 


’“1973 Voyaaer Capsule System Constraints and Required Docu- 
ment.” Internal document. 
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environment model developed for each atmospheric 
model, the points at which blackout b<‘pns and ends 
will be c*alciilated for 4()0 MHz and 2.295 GHz. These 
fr(H|uencies are the proposed transmission frequencies for 
communication link.s between the capsuh* and a relay 
bus, and directly between the capsule and the earth, 
respectively.’ 


will be encount(‘red by the capsule as it desc'ends through 
the planetars' atmosphere. Tlie properties of these mtxlel 
atmosph«*»’es are given in Table 1.* 

In the following, let h denote altitude above the Mar* 
tian surface. Then, the pressure p{h) of each atmospheric 
model is given by (Ref. 1) 


R. /v'artian Modal Atmospheres 

The physical properties of the Martian atmosphere and 
the variation of these properties with altitude above the 
Martian surface are not accurately known (Ref. 2). Recent 
information about the planet has led to the proposal of 
various model atmospheres. Since it is iina'rtain which 
model most accurately describes the actual Martian at- 
mosphere, two model atmospheres, the VM-4 and V’M-8, 
are used to develop two different plasma environment 
models to describe the actual pla.sma environment that 


Table 1. Propertiat of two modal .ifmotpharat 





Model 




VM-4 

VM-8 

Swffoc* praiiur* 

P" 

mbar 

10.0 

3.0 



Ib/ft’ 

20.0 

10.4 

Surface dentily 

P- 

iQ/cc) 10* 

2.37 

1.32 



(slugs/ft*) 10* 

4.98 

2.36 

Surface temperature 

r. 

•K 

200.0 

200 0 



•« 

360.0 

360.0 

Stratatpheric 

T. 

•K 

100.0 

100.0 

temperature 


•« 

180.0 

180.0 

Acceleration of gravity 

9> 

cm/s’ 

375.0 

373.0 

at lurfoce 


ft/s’ 

12.3 

12.3 

Compotition (percent) 





COi (by mou) 



70.0 

100.0 

COi (by volume) 



68.0 

100.0 

N] (by mou) 



0.0 

0.0 

N} (by volume) 



0.0 

0.0 

A (by moil) 



30.0 

0.0 

A (by volume) 



32.0 

0.0 

AAoleculor weight 

M 

mal'* 

42.7 

44.0 

Specific heat of mieture 

Cr 

cal/g*C 

0.1530 

0.166 

Specific heot ratio 

y 


1.43 

1.37 

Adiabatic lopte rote 

r 

*K/km 

-3.83 

-3.39 



•»/kft 

-3.21 

-2.96 

Trepopoute altitude 

h, 

km 

17.1 

18.6 



kft 

36.1 

61.0 


P 

km’’ 

0.193 

0.199 

(ttratoiphere) 


ft ' X 10* 

3.89 

6.07 


7 

ft/s 

133.3 

220.0 

wind tpeed 





Maximum surface 

V.u 

ft/s 

390.0 

336.0 

wind speed 





Design vertical 

dv/dh 

ft/s/kft 

2.0 

2.0 

wind grodient 





Design gust speed 


ft/s 

130.0 

200.0 


( r \T/«Y-n 

1 4- — /i j (1) 

for 0 ^ /i ^ h,, where p„ is the surface pressure, T„ is the 
surface temperature, y is the spec'ific heat ratio, r is the 
adiabatic lapse rate, and h, is the tropopause altitude, as 
given in Table 1. 

Similarly, the density’ p{h) and the surface density p„ 
of each atmospheric model are related by (Ref. 1) 

( p 

1 + /i j (2) 

for O^/i ^/i|. 

Above the tropopause altitude hi, each atmospheric 
model is assumed to be isothermal at the strato.«^pheric 
temperature of lOO’K. As a c'onsequence of this assump- 
tion, the pressure and density of each atmospheric model 
obey an exponential decay law with altitude of the form 
(Ref. 3) 

p = p.e/»* (3) 

P = P» ^ (4) 

for h> hi. The terms p, and pt art defined by 



(5) 

Pi = pi 

(6) 


where p< and p, are the pressure and density at the tropo- 
pause altitude as given by Eqs. (1) and (2), respectively. 

Figures 1 and 2 are graphs of pressure vs altitude and 
density vs altitude for the VM-4 and VM-8 model at- 
mospheres obtained from these equations. 
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Fig. 1. ^tme«ph*ric prctsur* vt oltitud* for 
two modoi otmosphorof 
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h, km 


Fig. 2. Atmotphoric dontity vt altitudo for 
two modoi atmotphorot 


III. Entry-Trajectory Characteristics 

A. Initial Conditions 


pt‘n<!mt of the* clisplac'ment of the capsule from the 
asstimtHl point of entry into atnv'rphere (Fig. 3). Clearly, 
the accuracy of this linear app.oximation increases with 
increasing entry angle. 


With tliese simplifying assumptions, the velocity v(h) 
and the elapst'd time-from-entr>’ it (h) of the capsule are 
given resp<*ctively by (Ref. 6) 

V ~ v,e ^ (7) 


» = — ^ r K( („ ) - Fi (, p: -)] 

^c,sin^, L \2g.Asin^,/ \2g.Asin^,/J 


(8) 


where p, is the pressure at the entry altihide, A is the 
ballistic c*oefficient of the capsule, and Ei (C) is the expo- 
nential integral function defined by (Ref. 5) 


for C > 0* 


EiU) 





£! 

V 


Figures 4 and 5 are graphs of velocity vs altitude and 
time-from -entry vs altitude for the \’M-4 and VM-8 mode) 
atmospheres obtain(*d from E<is. (7) and (8) for an entry 
velocity of 22 kft/s, an entry angle of 55 deg. and a 
ballistic coefficient of 0.12 slugs ft^ 


The entry phase of a Mars mission will begin at an alti- 
tude of 800 kft above the Martian surface. A set of initial 
conditions consistent with spacecraft approach trajec- 
tories. deflection maneuver orientation, and deflection 
maneuver accuracies is used for the entry-trajectory analy- 
sis. These conditions* at the entry altitude h, of 800 kft 
are (1) entry velocity o, of 22 kft s, and (2) entry angle 
of 55 rt6 deg. 

B. Entry-Trajectory Analysis 

The approach to the entr> -trajectory analysis is to use 
linear entry-trajectory theory to approximately determine 
the entry characteristics of the capsule. The entry char- 
acteristics of interest in this report arc the velocity and 
elapsed time-from-entry of the capsule calculated as a 
function of altitude. 

The entry-trajectory is assumed to be a straight path at 
the constant entry angle relative to the local horizontal 
at the assumed point of entry into the atmosphere. The 
surface of Mars is assumed to be flat so that the altitude 
of the capsule above the surface of the planet is inde- 


An examination of Fig. 4 reveals that, for the VM-4 
model atmosphere, the capsule achieves terminal velocity 



VERTICAL 


Fig. 3. Capsule entry-trajectoiy 
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Fig. 4. Capsule velccity vs altitude for two model atmospheres 


at approximately 50 s after entry. When this condition 
exists, h'(|. (8) is no longer valid and results in the iinrea- 
sonahly high value of 755.0 s to impact. Fortunately, only 
that part of the graph that lies below 50 s is of interest 
in the work that follows. 

coinputer program fo; the calculation of linear entr> - 
trajectories is described fully in tbe Append.x. 

IV. Plasma Interaction 

The interaction of an electromagnetic wave with a 
plasma can be di'scribed in terms of the parameters of 
the plasma. The parameters of interest in the present case 
are the electron concentration and the c'ollision frequency 
of the plasma. The.se parameters in turn depend on the 
interactions of the basic particles constituting the plasma. 
The number of free electrons contained within the plasma, 
or the electron concentration, is determined by the degree 
of ionization of the gases constituting the plasma. The 
number of collisions an electron undergoes per unit time, 
or the collision frequency, is determined by the inter- 
action of the free electrons with iens, neutral particles, 
and other electrons. 

From a knowledge of the spatial and temporal varia- 
tion of these parameters and a suitably chosen model 


to repre.sent the plasma, tbe frequency, permittivity, per- 
meability, and conductivity of the plasma can be deter- 
mi.ied. Tbe permittivity, }^*rmeability, and conductivity 
describe tbe macro.scopic behavior of the plasma and are 
called the constitutive parameters of the plasma. 

The interaction of an elc^ctromagnetic wave with the 
plasma can be determined from tlie propagation constant 
of the plasma. The propagao’on constant is a function of 
the plasma frt'tpiency and the constitutive parameters of 
the plasma and can be determined eithc*r by inverting the 
dispersion relation or by examining the wave eejuation in 
normal form. The real and imaginary parts of the prop- 
agation constant contain information alK)ut the phase 
shift and attenuation of the wave as it propagates in 
the plasma. 

The propagation constant of * plasma changes con- 
sideraldy depending on the relationship of the fr(‘tjuc*ncy 
of the electromagnetic wave in the plasma to the be- 
(iuency t>f the plasma. For frecjuencies above tbe plasma 
frequency, the plasma behaves as a lossy dielectric, and 
the electromagnetic wave, while experiencing some atten- 
uation, propagates freely. At frequencies below the 
plasma frecjuency, the plasma behaves .is a conductor, 
and the electromagnetic wave is highly a tenuated. 
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Fig. 5. Capsule fime-from-entry vs altitude for two model atmospheres 


Therv'fore, the interaction of an electromagnetic wave 
with the plasma can he determined, in principle, from a 
study of the relationships between the propagation con- 
stant, the constitutive parameters, the plasma freqiieiicy, 
the collision frecpiency, and the electron concentration of 
the plasma. In the following sections, these relationships 
arc developed; and, in particular, the propagation con- 
stant of the plasma is determined as a function of electron 
concentration and collision frecjuency. 

A. Constitutive Purameters of a Plasma 

suitable model of the plasma which is consistent with 
the limitations of this study is that of a certain number v 
of electrons per unit volume free to movs in an applied 
electromagnetic field, but subject to a damping force due 
to collisions characterized by damping conitant u>c. The 
flamping constant wc represents the average number of 
collisions the electrons undergo pe*^ anit time. The macro- 
scopic equation of motion of such electrons is 


nm ^ = nq (^' I v A;^) — nnu„,\ (9) 

where the applied electromagnedc field is character- 
ized by the field vectors ^ (r,t) and In the pres- 

ent case, the nonlinear v A •'JS term is dropped since 
|v A,:^^| << |<?|. 

In the steady state, with e **"' time-dependence, the 
equation of motion (9) becomes 

— io>nmv = nqE — nm<,>c\ (10) 

where the applied electromagnetic field in the steady state 
is characterized by the single field vector E (r). 

A rigorous derivation of the equation of motion using 
a statistical distribution function to describe the state of 
the plasma can be found in Ref. 6. 
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If the cximplex electric current density j(r) and 
radian plasma frequency are defined by 

the 

Let E (r) and II (r) be the complex electric and mag- 
netic field vectors in the plasma defined respectively by 

j M nqv 

(11) 

E(r)-He{<f (r,7)e 

(18) 

m€„ 

(12) 

H(r) — Re {,>/'(r,t)e 

(19) 


th(*n with the solution of E)q. (10) for v und the substitu- 
tion of the result into Eq. (11), the following relationship 
between j (r) and E (r) can be found: 


i(r) = 


r- / X / €,,o*,Wp 

n: — i 

loi T uir ^ ***? 



Tlie forn* of Eq. (13) suggests that the plasma has a 
complex conductivity ac given by 


where <f(r, f) and,X<'(r, t) are the real electric and mag- 
netic field vt*ctors in the plasma, respectively. For a linear, 


isotropic, and homogeneous plasma. 


D(r) = €E(r) 

(20) 

B(r) = MoH(r) 

(21) 

j(r) = aE(r) 

(22) 


For the avoidance of complex condtf;tivities, the plasma 
is considered to he a lossy dielectric. In this case the con- 
stitutive parameters are real, and they are found by sub- 
stituting Eq. (13) for j(r, /) into Maxwell’s equations. If 
the resulting equations are cast into the standard form 
for a lossy dielectric, it follows that the conductivity of 
the plasma ii given by 


where the complex vectors D (r), B (r), and j (r) represent 
the electric flux density, the magnetic flux density, and 
the electric current density of the plasma, respectively, 
and €, ft, and a are the constitutive parameters of the 
plasma as found in the previous section. 

With e time-dependence, .Maxwell’s two independ- 
ent equations take the following steady state form in 
the plasma: 

VA E(r) = i.,,/x„H(r) (23) 


a 


<!>■ + W? 


(15) 


VA H (r) = (a -«.,.€) E(r) (24) 


its permittivity is given by 



and its permeability is given by 




B. Propagation Constant of a Plasma 


(16) 


If Eq. (24) is substituted into the curl of Eq. (23) and 
the resulting curl-curl operation is expanded, the wave 
ecpiation for E (r) can be written in the form 

(V*+ y^’)K(r)=0 (2.5) 


(17) 


where the propagation constant y is given by 



(26) 


Let a plasma be characterized by its electron concen- 
tration n (r, t), its collision frequency ft- (r, t), and its per- 
meability ft = fto. 

If n(r, f) and /c(r, 0 are slowly varying functions of 
position in the interior of the plasma, then the field (juan- 
tities are continuous and have continuous derivatives and, 
therefore, satisfy Maxwell’s equations. 


The eonstant y is complex and, accordingly, is written 
in the form 

y = /34 ia (/?, a = positive-definite) (27) 

which displays the phase factor jS and the attenuation 
factor a. Explicit expressions for /? and a in terms of the 
eonstitutive parameters are found by substituting Eq. (27) 
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into Eq. (26). It follows that 
and 

To simplify the following discussion, it is assumed that 
the plasma is lossless. Since the losses are small (Kef. 3), 
this assumption is reasonable and intr(xluc(‘S negligible 
error. For a lossless plasma with a = 0, the expnvssions 
for P and « l)ecome 

= ' l‘l)l (30) 

o = “ + l‘l) J (31) 

Since € = €„(! — wJ/o>*) for a lossless plasma, 

P — ~ * TT = 0 for Wp < W 

(32) 

and 


/3 = 0 « = 


where the speed of light c in vacuum is defined by 


1 

V (mo^o) 


(34) 


Therefore, if wp < w, the wave travels without attenua- 
tion, and the plasma is said to be underdense. However, 
if u)p > <ij, the wave is evanescent and carries no power, 
and the plasma is said to be ovcrdensc. At o>p = w, the 
properties of the plasma change from underdense to over- 
dense, a!id the wave is cut off. This frequency is called 
the critical radian plasma frequency u»pc. 



f, Hz 


Fig. 6. Critical electron concentration vs frequency 


critical radian plasma frequency wpc. The critical electron 
concentration n, in terms of the signal frequency / is 
therefore 



Since the radian plasma frequency and the electron Loss of the transmitted signal will occur if the electron 
cxincentration are related by Eq. (12), a critical electron concentration exceeds the critical electron concentration 

concentration ric can be found that corresponds to each during entry into the Martian atmosphere. 
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Fig. 7. Plasma environment model 


From Kq. (35), the critical electron concentration for 
4(K) MHz is 1.99 X 10’ e*/cc, and tlie critical electron con- 
centration for 2.295 GHz is 6.53 X 10'”e*/cc. The values 
of the critical electron concentration for frequencies other 
than 100 MHz and 2.295 GHz are shown graphically in 
Fig. 6. 

C. Electron Concentration 
/. Stafituition region 

.An e.xact analysis of the electron conce ntration in the 
stagnation region (Fig. 7) of the capsule during (*ntry into 
the Martian atmosphere is not attempted here because 
the analysis would require the use of nonequilibrium 
chemistry to obtain the .solution. Instead, this analysis 
uses an appro.\irnate solution based on chemical equilib- 
rium in the stagnation region and is applicable only in 
the dense region of the atmosphere where the chemical 
reaction rates are sufficiently fast for the assumption of 


chemical eipiilihrium to be reasonable (Ref. 3). However, 
an equilibrium solution at less dense altitudes would pro- 
duce a conservative estimate oi the electron concentration 
and, therefore, will be in keeping with the objective of 
this study. 

With the use of an ( quilibrium thermochemistry and 
normal-shock computer program (Ref. 7), the equilibrium 
electron concentration behind a traveling normal shock is 
determined. This is a close approximation to the actual 
electron concentration in the stagnation region of a blunt 
entry capsule. 

The input parameters to the program are 
initial mixture constituents 
initial mole fraction compositions 
initial molal enthalpies 
initial free stream temperature 
initial free stream pressure 
shock temperature. 
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Tlu* output of tiu* program is a complftt* tluTmo<ly- 
numic aiui cliomical description of the gases behind the 
normal slunk. 

Tlie output parameters of interest to this report are 
the sho< k velm ity and electron amcentration behind the 
normal slunk. Since it is convenient to plot the Oiiti)iit 
parameters against altitude and since the program input 
is expressed in terms of initial free stream pressure, 
E(js. (1) and (3) are used to convert each altitude of 
interest into an ecjuivalcnt initial free stream pressure. 

Since the input value of the slunk temperature cannot 
he determined in advance, the electron c^onerntration in 
the stagnation region of the capsule is calculated for a 
series of shock temiXTatures. The program output then 
giv(‘S the appropriate electron concentration and shock 
velocity corresponding to each shock temperature in the 
series, from which a curve of electron concentration vs 
shoek velocitv' is developed for each initial free stream 
pressure, with the shock temperature as a parameter 
along the curve. The correct shock velocity for each curse 
is the velocity of the capsule for that curve; and the cap- 
sule velocity corresponding to the free stream pressure 
along each curve is given by Eq. (7). Equation (7) can, in 
turn, be used for the determination of the electron concen- 
tration in the stagnation region of the entr>’ capsule as a 
function of capsule velocity by linear interpolation, when 
necessary. 

Figure 8 shows graphs the electron concentration in 
the stagnation region of the capsule vs shock velocity or 
capsule velocity. These graphs were obtained for the 
V^M-4 and VM-8 model atmospheres by the method just 
described from tbe thermochemistry and normal-shock 
computer program. The cur»’es used to construct these 
graphs are also shown. 

2. Wake region 

The electron concentration in the wake region (Fig. 7) 
is calculated since the antenna will be mounted on the 
rear part of the capsule and the propagation path will be 
througb the wake region (Ref. 1). 

.\gain, no exact analysis is attempted and a frozen flow 
approximation to the actual flow conditions is used to 
estimate the value of the electron concentration in the 
wake region of the capsule (Ref. 3). It is assumed that the 
total number of electrons formed in the stagnation region 
remains c-onstant as the gases flow around the capsule and 
into tbe wake region. .Any change in the electron concen- 
tration is due to the expansion of the gases as they flow 


into the wake region. Tlie value of the specific heat ratio 
is not a constant during the expansion of the gases; how- 
ever, the results are not sensitive to the value selected 
(Ref. 3). 

The value of the electron exmeentration in the wake 
region is related to the value of the eh'ctron concentration 
n, in the stagnation region by (Ref. 5) 



where the pressure ratio p,/p« across the sh(K-k is ob- 
tained from the output of the thermochemistry and 
normal-shoek computer program, and y is the specific 
heat ratio. 

Graphs of electron concentration in the wake region vs 
altitude, and electron concentration in the wake region 
vs time-from-entry, can be obtained from Fig. 8 by using 
E(|. (36) and the relationships between capsule velocity, 
altitude, and capsule time-from-entry as given by Eqs. (7) 
and (8). These graphs are shoum in Figs. 9 and 10 for the 
\'M-4 and VM-8 model atmospheres. Note that the values 
of the electron concentrauon between 22 km and 100 km 
for the V^M-4 model atmosphere and bc*tween 20 km and 
80 km for the VM-8 model atmosphere were obtained 
from the thermochemistry and normal-shock computer 
program. 

V. Conclusion 

.As is evident from an examination of the graphs of 
electron concentration in the wake region vs capsule alti- 
tude or time-from-entry, blackout will occur during entry 
of a blunt capsule into the Martian atmosphere. 

The electron concentration in the wake region reaches a 
peak of 7.7 X 1 1 e*/cc at an altitude of 29 km for the VM-4 
model atmosphere and reaches a peak of 7.9 X 10’® e'/cc 
at an altitude of 29 km for the VM-8 model atmosphere. 

In the VM-4 model atmosphere, blackout begins at an 
altitude of 66 km and ends at an altitude of 22 km for the 
400-MlIz signal, and begins at an altitude of 47 km and 
ends at an altitude of 24 km for the 2.295-GHz signal. The 
duration of blackout is 8.5 s for the 400-MHz signal and 
4.5 s for the 2.295-GHz signal. Similarly, in the VM-8 
model atmosphere, blackout begins at an altitude of 
53 km and ends at an altitude of 23 km for the 400-MHz 
signal, and begins at an altitude of 33 km and ends at an 
altitude of 27 km for the 2.295-GHz signal. The duration 
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of blackout is 5.7 s for the 4(X)-\IHz signal and 1.3 s for 
tlu* 2.295-CH/. signal. 

Blackout can he rcdiic'cd or climinatt'd by entering the 
atmosphere at lower velocities, by transmitting at higher 


frequencies, or by injecting high-f'lectmn-affinity fluids 
into the wake region. TIk* latter alternative has b<*en 
examined by D. F. Spencer (Ref. 3} anti the feasibility of 
the first and seconti alternatives will be the subjt'cts of 
future invi*stigations. 
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Appendix 

Th« Calculation of Linear Entry-Trajoctorlos: 
A Computor Program 


I. Progrcim Description 

The purposf of this pro^am, which has liecn speci- 
fically developed for use with the Martian entry propa* 
gation study, is to determine the linear entry-trajectory 
characteristics of a blunt capsule as it dv^scends through 
a Mars-type planetar)’ atmosphere subject to arbitrary 
initial entr>' amditions. Acttially, the development of 
the program is such that any hypothetical planetary 
atmosphere can be use<! in the calculations of the entry- 
trajectory characteristic's of the capsule. Tfie character- 
istics of interest to this study are the veltK-ity and elapsed 
time-from-entry of the capsule as a function of attitude. 
The source language is FORTRAN II for use on an IBM 
1620 computer or equivalent. 

A. Input Raquiramants 

To compute the entry-trajc*ctory characteristics, the 
pressure and density profiles of the atmosphere must first 
be developed. This program develops these profiles from 
the assumt'd physical properties of the atmosphere. The 
requirt'd input data to the program are 

(1) Atmospheric properties 

surface pressure 
surface density 
surfac'e temperature 
specific heat ratio 
adiabatic lapse rate 
tropopause attitude 
inverse scale height 

(2) Entry parameters 

entry attitude 
entry velocity 
entry angle 
ballistic c'cx'fficient 


B. Output 

The output of the program includes the pressure and 
density profiles of ihe atmosphere and the velocity and 
*ime-from-entry of the capsule. The various outputs of 
the program are calculated at points from the surface of 
the planet to any arbitrary altitude above the surface 
to the planet in arbitrary increments in altitude. 


C. Program Limitotions 

Although the program has certain limitations, they are 
entirely in keeping with the study. 

The limitations follow. 

(1) TIu‘ trajc*ctory is assumed to be a straight path 
at the constant entry angle relative to the local 
horizontal at the assumed point of entry into the 
atmosphere. 

(2) The surface of the planc't is assumed to be flat so 
that the altitude of the cap.sule above the surface 
of the? planet is indepc'ndent of the displacement 
of the capsule from the assumt'd point of entry into 
the atmosphere. 

(3) Above the tropopause altitude, the planetary atmo- 
sphere is assumt'd to be isothermal at the strato- 
spheric temperature so that the pressure and density 
profiles obey simple exponc'iitial det»y laws. 


II. Calculations 

A. Atmospheric Charactoristics 

The pressure and density profiles of the atmosphere 
are given respectively by 

( r 

1 + y /i j O^h^ht 

h>ht 

and by 

1 / r 

h > ht 

where p„ is the surface pressure, p„ is the surface density. 
To is the surface temperature, y is the specific heat ratio, 
r is the adiabatic lapse rate, ht is the tropopause altitude, 
Pi is the pressure of the tropopause altitude, and pt is the 
density at the tropopause altitude. 
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ft. Er>try-Tro|«ctory Choroct«rittic« whrrr 

The* velocity i: and time-from-entry 8t of the capsule arr 
Ipvm respectively by 



lor C > 0 




and 


/^v, sin ( 2^',:^ sin ( 2fi,^ sin )] 


In these equations, p, is the pressure at the entry altitude, 
V, is the entry velocity, is the entry angle, A is the ballis- 
tic coefficient, R, is the surface gravity, fi is the inverse 
scale height, and Ei{C) is the ex{)onential integral func- 
tion. 


III. Usage 

A. Definitions 

Program usage is given in the form of tables that fol- 
low. Input and output data are shown in Tables A-1, A-2, 
an<l A -3. 


Tablo A-1. Input definitions 


Pr» 9 ram noma 

SymO*l 

Pr«9«rty 

Unit 

PO 

e- 

Swrfoc* pratiura 

mbor 

DC 

P' 

Surface daniily 

10 *g/cc 

TO 

T. 

Surface (emperaiure 

•K 

SC 

»• 

Surface gravity 

cm/i’ 

SHR 

y 

Specific heal ralia 

— 

AIR 

r 

Adiabatic lapie rate 

*K/lim 

MT 


Tropapaute altitude 

km 

SHI 



km’* 

HE 

h. 

Entry altitude 

kft 

VE 

¥, 

Entry velacily 

kfl/i 

EA 


Entry angle 

deg 

•C 

A 

Salliitic coefficient 

ilugi/ft’ 

AIT 

O 

Altitude to which calculatiani 
eitend 

km 

DELTA 

S 

Increment in altitude between 
calculatiani 

km 


Table A-2. Input data 


Cord number 

Program name 

Permot 

1 

PO. DO. TO. SO. SHR. AIR, HT, SHI 

• PIO.O 

2 

HE. VE, EA, K. ALT. DELTA 

6M0.0 


Table A-3. Output definitions 


Program name 

Symbol 

Property 

Unit 

H 

h 

Altitude 

km 

HC 

h. 

Altitude 

kft 

P 

P 


mbor 

PC 

Pf 

Preitura 

P/*’ 

PA 

P* 

Preiiure 

olm 

D 

P 

Dentily 

10 *g/cc 

DC 

P' 

Deniity 

10 *ilugi/f’ 

V 

V 

Velocity 

km/i 

VC 

Vr 

Velacily 

kfl/i 

T 

T 

Time-from-entry 

ft 
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B. Sample Output Data 

/. The VM-4 model a/mo».*/iere 


SURFACE PRESSURE 
Surface density 

SURFACE TEMPERATURE 
SURFACt gravity 
SPECIFIC HEAT RATIO 
ADIABATiC LAPSE RATE 
TROPOPAUSE altitude 
<NVERSE scale height 

ENTRY altitude 
ENTRY velocity 

entry angle 
ballistic coefficient 


10.000 

(M8I 

2.570 

IE-05 G/CCI 

200.000 

lOEG.K) 

375.000 

ICK/SO.SI 

1.430 
-5. *50 

IDEC.k/kmi 

17. 100 

IRM) 

.193 

1 1/KMI 

800.000 

IKF) 

22.000 

UF/SI 

55.000 

(DEGREES) 

.120 

1 SLUGS/ S';. F) 


altitude 

ALTI TUDE 

PRESSURE 

PRESSURE 

PRESSURE 

DENSITY 

density 

IKM) 

(KF ) 

(MB) 

(P/SO.F) 

(ATMSPHRS) 

IE-05 G/CC) 

(F-05 SLUGS/CU.E 


0.000 

0.000 

1 .OUOOOOUE^OI 

2.0685514E«01 

9.8692G00E-03 

2.5700000E>00 

4.98663086*00 

1.000 

3.280 

9.0599252E»00 

1 .89221 19E401 

8.9414213E- 73 

2.39855856*00 

4.65397806*00 

r.ooo 

6.561 

8. 1834675E^00 

1.7091592E^01 

8.0764277E-03 

2.‘2338302E*00 

4.33435276*00 

3.000 

9.842 

7.3680988E^U0 

1 .53HH653E«01 

7.2717240E-03 

2.07574836*00 

4.02762276*00 

4.000 

13.123 

6.611 3163E-»00 

1 .380H073E^01 

6.5248402E-03 

1 .92424496*00 

3.73365726*00 

5.000 

16.404 

5.9106429E«00 

1 .23446ilE^0l 

5.8333316E-03 

1.77925066*00 

3. 452^21 36*00 

6.000 

19.6B4 

5.2636285E>00 

1 .0993358E^0l 

5. 1947802F-03 

1.6406943E *00 

3.18347736»00 

7.000 

22.965 

4.6678501E^00 

9.7490448E^0U 

4.6067946E-0? 

1 ,50850356*00 

2.92698446*00 

8.000 

26.246 

4. 1209123E«00 

8.6067371E«0U 

4.0670107E-03 

1.38260376*00 

2.68269816*00 

9.000 

29.527 

3.6204484E-»00 

7.5614925E^00 

3.5730929E-03 

1 .2629185F*00 

2.45047016*00 

10.000 

32. BOH 

3. 1641213E'»00 

6.60H4299E«0U 

3. 1227\45E-03 

1 . 14936976*00 

2.23014886*00 

11.000 

36.089 

2. 7496226E^U0 

5.7427281E>0U 

2.7136575E-03 

1 .04187686*00 

2.02157786*00 

12.000 

39.369 

2.3746766E-»00 

4.9596341E^00 

2.3436158E-03 

9.40357346-01 

1.82459726*00 

13.000 

42.650 

2.0370384E-»00 

H .2544594E>0U 

2.0103939E-03 

8.44725876-01 

1.63904126*00 

14.000 

45.931 

1.7344956E^00 

3.6225832E-^0U 

1 .71 18083E-03 

7.54894766-01 

1 .46473986*00 

15.000 

49.212 

;.46<.8699E«’00 

3.0594560E^0U 

1.4457094E-03 

6.7077344E-01 

1 .30151736*00 

16.000 

52.493 

1.2260180E+00 

2.5606016E«^0U 

1 .2099816E-03 

5.92268146-01 

1.14919186*00 

17.000 

55.774 

1.0158325E^OO 

2 .12161 8^: »0U 

1 .0025454E-03 

5.19281816-01 

1.00757456*00 

18.000 

59.054 

8.3744778E-01 

1 .749052 7E *00 

8.2649396E-04 

4. 3059H89F-01 

8.35501056-01 

19.000 

62.335 

6.9046060E-U1 

1 .44201. 24E«^0U 

6.8142937E-04 

3.55021026-01 

6 . 888555S6-() 1 

20.000 

65.616 

5.6927231E-01 

1 .1889544E^0U 

5 .61H2622E-04 

2.9270H43E-01 

5.67948986-01 

21.000 

68.897 

4.6935474E-01 

9.8027149E-01 

4.6321558E-04 

2.41332826-01 

4.6826^696-01 

22.000 

72. 178 

3.8697454E-01 

8.U321621E-01 

3.8191 291E-04 

1 .98974566-01 

3.86074976-01 

23.000 

75.459 

3. 1905353E-U1 

6.6635969E-01 

3. 1468030E-04 

1 .64050946-01 

3.18311866-01 

24.000 

78.739 

2.6305389E-0I 

5.4940157E-01 

2.5961314E-04 

1.35257046-01 

2.62442386-01 

25.000 

82.020 

2. 1688319E-01 

4.5297169E-01 

2. 1404635 E-04 

1. 11516996-01 

2. 16379016-01 

26.000 

85.301 

1.7881628E-01 

3.7346699E-0I 

1 .7647736E-04 

9. 19437536-02 

1.78400606-01 

27.000 

88.582 

1.4743080E-01 

3 .U791680E-01 

1.4560240E-04 

7. 58059686-02 

1 .47088086-01 

28.000 

91.863 

i.2155403E-01 

2.5387183E-01 

1. 1996410E-04 

6.25006536-02 

1.21271476-01 

P'^.OOO 

95 . 1 44 

1 .0021911E-01 

2.O931276E-01 

9.8O08244E-05 

5. 1 53C6626-0? 

9.99861446-0? 

;n.000 

98,424 

8.2628H50E-02 

1 .7257460E-01 

8.1S48064E-05 

4.24861016-02 

8. 24367716-02 

;i.ooo 

101.706 

6.81259f'7E-02 

1 . 42284 64E-01 

6.7234908E-05 

3.50290246-02 

6.79676316-02 

32.000 

104.986 

5.6168657E-02 

1.1731112E-01 

5.5433970E-05 

2.88807996-02 

5.60380876-02 

33.000 

108.267 

4.6310045E-U2 

9.6720909E-02 

4.5704309E-05 

2.38116986-02 

4.62023926-0? 

34.000 

111.548 

3.8181796E-02 

7.9744643E-02 

3.7682378E-0? 

1.96323156-02 

3.80930386-02 

35.000 

114.829 

3. 1480203E-02 

6.5748022E-02 

3. in6844lf-05 

1 .61864906-02 

3. 14070236-02 

36.000 

118.109 

2.5954861E-02 

5.4208061E-02 

2.5615371E-05 

1 .33454706-02 

2.58945266-02 

37.000 

121.390 

2. 1399314E-02 

4 .4693567E-02 

2. 1119410E-05 

1 . 10030996-02 

2. 13495696-02 

38.000 

124.671 

1.7643349E-02 

3.6849041E-02 

1 .7412573E-05 

9.07185726-03 

1 .76023356-0? 

39.000 

127.952 

1.4546623E-02 

3.0381369E-02 

1.4366353F-05 

7.47958276-03 

1 .45128086-02 

40.000 

131.233 

1 . 1993428E-02 

2.5048890E-0^ 

1. 183655 3E -05 

6.16678126-03 

1 .19655496-02 

41.000 

134.514 

9.8883659E-03 

2 .0o5:360E-02 

9.7590260E-06 

5 .08439996-03 

9.86537956-03 

42.000 

137.794 

8. 1527T94E-03 

1 .70274V8E-02 

8.0461410E-06 

4. 19199606-03 

8. 13382756-03 

43.000 

141 .075 

6.7218196E-03 

1 .4038865E-0^ 

6.6338981E-06 

3.45622516-03 

6.70619406-03 

44.000 

144.356 

5.5420191E-03 

1.1574791E-02 

5 .4695294E-06 

2.84959536-03 

5.52913616-03 

45.000 

147.637 

4.5692949^-03 

9.5432072E-03 

4.50952856-06 

2.34944006-03 

4.55667316-03 

46.000 

150.918 

3.7673Oi2E-03 

7.8682021E-03 

3. 71802496-06 

1.93707096-03 

3.75854386-03 

47.000 

154.199 

3.1060718E-03 

6.4871906E-03 

3.06544436-' 6 

1 .59707996-03 

3.09885136-03 

48.000 

157.479 

2.5609003E-03 

5.3485719E-03 

2.52740376-06 

i.?l676376-03 

2.55494726-03 

49.000 

160.760 

2.1114162E-03 

4.4098012E-03 

2.08379886-06 

1.08564796-03 

2.10650706-03 

50.00'J 

164.041 

1.7408245E-03 

3.6358014E-03 

1 .7180545F-06 

8.95097186-04 

1.73677796-03 
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At T ITUOe 

AlTITUOf 

VFLtXITV 

VELOC ITV 

T IMF FAOM FMBV 


IKF 1 

(KM/SI 

<KF/S > 

( SFC. 1 


0.000 

0.000 

1. 1922950E-03 

3.91 1 721 IE 03 

7. *4981758+02 

1.000 

3.280 

2.6848309E-03 

8.H0H4826E-03 

3.98662358+02 

2.000 

6.561 

5.7226001E-03 

1.8774896E-02 

r .31596008+02 

3.000 

9.f!42 

1. 1570739E-02 

3.7961665E-02 

1.4882312E+02 

4.000 

13. 123 

2.2241225E-02 

7.2969751E-0C 

I .05580848+02 

5.000 

16.404 

4.0730062E-02 

1 .3362854E-01 

8. 18*0625E+01 

6.000 

19.684 

7.1211314E-02 

2.3363245E-01 

6.8040523E+01 

7.000 

22.965 

1. 1911574E-01 

3.9079888E-01 

5.97921458+01 

8.000 

26.246 

1.9101780E-01 

6.2669755E-01 

5.4543295E+01 

9.000 

29.527 

2.9427307E-01 

9.6546088E-01 

5. 10782288+01 

10.000 

32.808 

4.3639092E-01 

1.4317258E40U 

4.86953898+01 

11.000 

36.089 

6.241 8741E-01 

2 .0478548E^00 

4.69917158+01 

12.000 

39.369 

8.6282314E-01 

2.b307788E40U 

4.57276158+01 

13.000 

42.650 

1.1548829E400 

3.7889782E40U 

4.47560028+01 

14.000 

45.931 

1.4996612E^00 

4.9201 384E400 

4. 39837458+01 

15.000 

49.212 

1.8927980E+00 

6.2099547E400 

4.37500998+01 

16.000 

52.493 

2.3263498E+C0 

7.6323658E400 

4.28142568+01 

17.000 

55.774 

2.7893033E400 

9.1512391E400 

4.23479648+01 

18.000 

j9.054 

3.2538011E400 

1 .0675179E*01 

4.19431268+01 

19.000 

62.335 

3.6941285E400 

1 .2119819E401 

4. 1591493E+01 

20.000 

65.616 

4. 1016469E+00 

1 .3456819E+01 

4.12782598+01 

21.000 

68.897 

4.4712392E400 

1 .4669390E401 

4.09935028+01 

22.000 

72. 178 

4.8008809E400 

1 .5750889E401 

4.07302608+01 

23.000 

75.459 

5.0908662E400 

1 .6702283E401 

4.04835258+01 

24.000 

78.739 

5.3430822E+00 

1 .7529762E401 

4.02496148+01 

25.000 

82.020 

5.5604011E400 

1 .b242749E+01 

4.0025V75E+01 

26.000 

85.301 

5. 7462091E400 

1 .bH«,2354E401 

3.9809910E+01 

27.000 

88.582 

5.9040652E>00 

1 .9370253E401 

3.9600405E-01 

28.000 

91.863 

6.0374713E400 

1 .9807936E401 

3.9396003E+01 

29.000 

95.144 

6. 1497263E+00 

2 .0176226E+01 

3.91957148+01 

30.000 

98.424 

6.2438467E-^00 

2.U4B5020E401 

3.8998753E+&1 

31.000 

101 .705 

6.3225298E400 

2 .0743166E401 

3.88044958+01 

32.000 

104.986 

6.3881478E400 

2.0958448E401 

3.86124378+01 

33.000 

108.267 

6.4427605E400 

2.1137623E401 

3.84221748+01 

34.000 

111.548 

0.4881386E4U0 

2.1286501E401 

3.82333788+01 

35.000 

114.829 

6.5257925E400 

2 .1410037E401 

3. 80457858+01 

36.000 

118.109 

6.5570016E+00 

2.1512429E401 

3.78591758+01 

37.000 

121.390 

6.5828452E400 

2 .1597217E401 

3.76733748+01 

38.000 

124.671 

6.6042294E400 

2.1667375E+01 

3.74882368+01 

39.000 

127.952 

6.6219126E400 

2 .1725391E+01 

3.73036448+01 

40.000 

131.233 

6.6365276E4C0 

2.1773340E+01 

3.71194998+01 

41.000 

134.514 

6. 648601 7E400 

2.1812953E+01 

3.69357248+01 

42.000 

137.794 

6.6585731E400 

2.184566BE+01 

3.6752252E+01 

43.000 

141 .075 

6.666805£E400 

2.1872677E+01 

3.65690308+01 

44.000 

144.356 

6.6736007E+00 

2.1894971E401 

3.63860148+01 

45.000 

147.637 

6.6792084E4.00 

2.191336VE401 

3.62031678+01 

46.000 

150.918 

6.6838'’ 4E+00 

2 .192854 9t +01 

3.60204608+01 

47.000 

154.199 

6.6876527C400 

2 .1941073E+01 

3.58378688+01 

48.000 

157.479 

6.6908017E400 

2.1951405E+01 

3.56553718+01 

49.000 

160.760 

6.693399CE400 

2.1959926E+01 

3.5472952E+01 

50.000 

164.041 

6.69554l2E-*00 

2.1966954E+01 

3.52905978+01 
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2. The VM-H model atmosphere 


SURFACE PRESSURE 
SURFACE OcNSITV 
SURFACE TEMPERATURE 

surface gravity 

SPECIFIC MEAT RATIO 
ADIABATIC LAPSE RATE 
TROPOPAUSE altitude 
INVERSE scale HEIGHT 

entry altitude 

ENTRY VELOCITY 
ENTRY ANGLE 
BALLISTIC COEFFICIENT 


5 

.000 

«MBI 

1 

.320 

(E-05 G/CCI 

?00 

.000 

(OEG.KI 

375 

.000 

ICM/SU.SI 

1 

.3 70 


-5 

.390 

(DEG.K/KM) 

IH 

.F>00 

IKMI 


.199 

I 1/KMI 

800 

.000 

(KF I 

22 

.000 

(KF/SI 

55 

.000 

(DEGREES) 


.120 

ISLUGS/SU.F) 


altitude 

altitude 

PRESSURE 

PRESSURE 

PRESSURE 

DENSITY 

DENS I TV 

(KM) 

(KF ) 

(MB) 

(P/SU.F ) 

( ATMSPhrs ) 

(E-05 G/CC) 

(F-05 SLItGS/CU.F ) 


0.000 

0.000 

5.0000000E^00 

1 .0662757E401 

6.9366000E-03 

1 . 3200000E+00 

2.5M 2267F + 00 

l.COO 

3.280 

6.51R9550E-*-00 

9.6380697E»0U 

6.659H670F-03 

1 .2260<*60F4’00 

2.37h92S6F<-00 

2.000 

6.561 

6.07?6025E^00 

8.5058396E+0U 

6.019332B*--03 

1. 1366?niF^00 

2.205(i22^k^C0 

3. DOC 

9.862 

3.6593226E400 

7.6626H33E^0y 

3.61I6586F-03 

1 .0510375F*00 

2.0393525F400 

A. 000 

13.123 

3.2775276E+U0 

6.H652 866E<'00 

3.2366573F-03 

9.69H1309F-01 

I .88 175098 *00 

5.000 

16.606 

2.92566l5E<-()0 

6.1 103966E*0u 

2.P873938E-03 

8 .9266065F-01 

1 .732ii5n?F<.()0 

6.000 

19.686 

2.6022020E+00 

5.636B326E-»0U 

2.56H1651E-03 

8. 1969362F-01 

1 .59u0821F^OO 

7.000 

22.965 

2.3056592E-»00 

6.B156877E+00 

2.2755011 E-03 

7.5022372F-01 

1 .65*>6765F>0(j 

B.OOO 

26.266 

?.0365766E^00 

6.269317hE*0O 

2.0079663E-03 

6.8676316E-01 

1.32R6618F*00 

9.000 

29,527 

1.7875313E+00 

3.73335C9E^OU 

1.7661503E-03 

6.2302229E-01 

1 .20rth666F4.00 

10.000 

32.B08 

1.5631362E+00 

3.26668tiC t-OU 

1 .5626886F-03 

5.66910B8F-01 

1 .096l097F + 0(» 

11.000 

36.089 

1.3600307E^OO 

2 .B606960E^0(J 

1.3622616E-03 

5. 103376MF-01 

9.9022O09F-O1 

12.000 

39. 369 

1. 1769007E+00 

2.6580176E<-0w 

1. 1615068E-03 

6.5921066E-01 

M.91U1671F-U1 

13.000 

62.650 

1 .0126595E+00 

2 .1 165737E+OU 

9.9921652F-06 

6.116358HF-01 

7.9R31R62F-01 

16.000 

65.931 

8.6565765E-01 

1 .H075527E^0U 

8.S613766E-06 

3.6691958E-01 

7.11962MF-C1 

15.000 

<•9.212 

7.366B160E-01 

1 .5366198E*0u 

7.2S07176E-06 

3.2556593E-01 

6.3l7t>3l5F-ni 

16.000 

52. < 3 

6. 18953B0E-01 

1 .292716BE+0U 

6. 1OH57H8E-06 

2.B727816F-01 

5.5761766F-01 

17.000 

55. , 76 

5. 1 713630E-01 

1 .08006 15E+00 

5.1 370188-06 

2.51O5806E-01 

6.PrtBH00SF-Ul 

IB. 000 

59.056 

6.2812006E-01 

H.9615075E-01 

6.2252026F-06 

2. 195061 lE-01 

6.2591281F-01 

19.000 

62.335 

3.5130627E-01 

7.3372120E-01 

3.6671 1 18E-06 

1 .8S96205E-01 

3.6062664P-01 

20.000 

65.616 

2.8791301E-01 

6.0132112E-0I 

2.861<.7inE-06 

1.5260518F-01 

2.957) 532F-01 

21.000 

68.897 

2.3595907E-01 

6.9281266E-01 

2.3?87272E-06 

1 .2690366E-1)1 

2.6235766F-01 

22.000 

72. 178 

1 .9338023E-U1 

6.0388655E-01 

1 .9085081E-06 

1.0236677E-01 

1 .98o207<.t -01 

23.000 

75.659 

1 .586B675E-01 

3.3100356E-01 

1.5661176F-06 

H.3R93039F-02 

1.62 /7961F-01 

26.000 

78.739 

1.2988616E-U1 

2.7127392E-01 

1.2818726F-06 

6. H756532F-02 

1 .3360602F-01 

25.000 

82.020 

1 .0666819E-01 

2 .2232251E-01 

1 .0505586F-06 

5.b367769E-02 

1 .U9332IIRF-01 

26.000 

85.301 

8.7239606E-02 

1 .8220639E-01 

8 .6098S09E-05 

6.6179H08F-02 

H.96U3767F-02 

27.000 

88.582 

7. 1697209E-02 

1 .6932559E-01 

7.0567025F-05 

3.7866657F-02 

7.363‘.760F-02 

2B.000 

91.863 

5.8595529E-02 

1.2237977E-01 

5 .7829099E-05 

3. lO17223fc-02 

6.01836605-02 

29.000 

95. 166 

6.8021960E-02 

i .0029633E-01 

6.7393832E-05 

2.5620162F-02 

6.U323332F-02 

30.000 

98.626 

3.9356392E-02 

8.2197H67E-0^ 

3.8B61610E-05 

2.0833091F-02 

6.062/9315-02 

31.000 

101.705 

3.2266526E-02 

6.7365235F-0^ 

3. 1832636E-05 

1.7073757F-02 

3 . 3128608F-(12 

32.000 

106.986 

2.6636193E-02 

5.5209l70E-0i 

2 .608P633F-05 

1.3992795E-C2 

2.71505665-02 

33.000 

108.267 

2. I666161E-02 

6.5266672E-0^ 

2. 138O776F-05 

1 . 1667796E-02 

2.2251227F-U2 

36.000 

111.568 

1 . 7756866E-02 

3.70819068-02 

1 .7522610E-05 

9. 3986299E-03 

1 .6235992F-U2 

35.000 

116.829 

1.6550981E-02 

3.0390671E-0^ 

1.6360656F-05 

7.7026H23E-03 

1 ,4965306F-02 

36.000 

118.109 

1. 1925256E-02 

2.6906505E-02 

1.1769271E-05 

6.3125685E-03 

1.2268623F-02 

37.000 

121.390 

9.7733616E-03 

2 .0612125E-02 

9.6655060E-06 

5.1736669F-03 

1.00351565-02 

3B.000 

126.671 

8.0097609E-03 

1 .6728755E-02 

7.9069736E-06 

6.2399126F-03 

8.22650125-03 

39.000 

127.952 

6.5663822E-03 

1 .3710069E-02 

6.6785200F-06 

3.676H19HF-03 

6.762273HF-03 

60.000 

131.233 

5.379B387E-03 

1.1236069E-02 

5.3096706E-06 

2.R6778H2F-03 

5.52562975-03 

61.000 

136.516 

6.6090662E-03 

9.2085196E-03 

6.3513758E-06 

2.3339065E-03 

6.52852925-03 

62.000 

137.796 

3.6136336E-03 

7.566861 3E-03 

3. 5661696k -06 

1.9127512E-03 

3i71 13557F-03 

63.000 

161.075 

2.9613889E-03 

6.1850129E-03 

2.9226539E-06 

1 .5675950E-03 

3.06166115-03 

66.000 

166. 356 

2.6270060E-03 

5.0689267E-03 

2.3952607E-06 

1.2867223E-03 

2.69277665-03 

65.000 

167.637 

1.9R90525E-03 

6 .1562383E-03 

1.9630356E-06 

1 .0528960E-03 

2.06295675-03 

66.000 

150.918 

1.6301277E-03 

3 .6066056E-03 

1 .6088056E-06 

8.6289923E-06 

1 .676?0365-( 3 

67.000 

156.199 

1 .3359710E-03 

2 .7902661E-03 

1. 3186966E-06 

7.0718898F-06 

1.37217525-03 

68.000 

157.679 

1 .0968969E-03 

2.2867662E-03 

1 .0805736E-06 

5.7O57665F-06 

1 .12656605-03 

•♦9.000 

160.760 

B.9732103E-06 

1 .8761010E-03 

8.8558607F-07 

6.7‘.9919HF-06 

9. 21638005-06 

50.000 

166.061 

7.3539936E-06 

1 .5359193E-03 

7,2578031E-07 

3. 8927962F-06 

7.55328315-06 
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LT ITUDE 

alti T unp 

vFLOC nv 

VELOC ITV 

TIME FROM ENTRY 

(KM» 

«KP ) 

(KM/S> 

(KF/S 1 

(SEC. ) 


0.000 

0.000 

H.9415150E-02 

2. 933562 OE -01 

6.3934973E^01 

1.000 

3.280 

1.3545829E-01 

4.4441606E-01 

5.8275553E>01 

2.000 

6.561 

1 .9915457E-01 

6.5339294E-01 

5.4380182E>01 

3.000 

9.842 

2.8455912E-01 

9.3359103E-01 

5. 1616908E^01 

A. 000 

13.123 

3.9568338E-01 

1.2981 712E*00 

4.o599415E^01 

5.000 

16.404 

5. 3i ’ 661 8E-01 

1 .7590718E^0U 

4.8085368E-»01 

6.000 

19.684 

7.0892207E-01 

2.3258551E^0U 

4.6918960E*01 

7.000 

22.965 

9. 158064 7E-01 

3 .U0460B3Et0U 

4.5997658E-*'01 

8.000 

26.246 

i. 1573404E-»00 

3.7970409E^00 

4.5252470E^01 

9.000 

29.527 

1.4325325E-»00 

4.6999003E*0U 

4.4635969E^01 

10.000 

32.808 

1. 738H200E+00 

5. 70477858^00 

4.4114883E^01 

1 1.000 

36.089 

2.0721422E+00 

6.7983531E>0o 

4. 366540 lE^Ol 

12.000 

39.360 

2.4271372F+00 

7.9630325E+00 

4.3270136E^01 

13.000 

42.650 

2. 7974391E-^00 

9.1779313E+0U 

4.2916148E+01 

lA.OOO 

45.931 

3. 1760486E+00 

1 .U4200«6E>01 

4.2593600E-f01 

15.000 

49.212 

3.555731 5E-»o0 

1 . 1665762E+01 

4.2294H50E^01 

16.000 

52.493 

3.9294108E*00 

1.2891 741E+01 

4.2013827F-f01 

17.000 

55.774 

4.2905237E*00 

1 .4076493E+01 

4. 17455 72E+01 

1 8 . 000 

59.054 

4. 6333076fc^00 

1.5201109E^01 

4. 1485920E+01 

19.000 

62.335 

4.9510434E+U0 

1 .6243548E-*’01 

4. 1232874Et.01 

20.000 

65.616 

5.22O6130E+O0 

1.7157488E>01 

4 .0993149F.f01 

21.000 

68.897 

5.4695621E-^U0 

1 .7944721E+01 

4.0768042E^01 

22.000 

72.178 

5.6743988E+00 

1 .8616756E^01 

4.0546035E+01 

23.000 

75.459 

5.8479797E>00 

1 .91H6246E+01 

4.0334213F-*'01 

24.000 

78.730 

5.9941 897E+00 

1 .9665937E+01 

4.0128101E+01 

25.000 

82.020 

6. 1 167376E-^U0 

2.0067996E+01 

3.9926556E+01 

26.000 

85.301 

6.2190375E+00 

2.0403625E-^01 

3.9728678E+01 

27.000 

88.582 

6.3041518E+00 

2.0682871E+01 

3.9533753E+01 

28.000 

91.863 

6.3747750E+U0 

2.0914574E>01 

3.9341217E>01 

29.000 

95.144 

6.4332437E+OU 

2.11064O0E-H. 1 

3.9150617E+01 

30.000 

98.424 

6.4B15614E+00 

2.1264922E+01 

3.3961587E+01 

31.000 

101.705 

6.5214307E-^00 

2.1395727E+01 

3.8773837E;-01 

32.000 

104,986 

6. 55428R2E+U0 

2.1503526E+01 

3.8587128E+01 

33.000 

108.267 

6. 5813401E+0U 

2.1592279E401 

3.8401268E^01 

34.000 

111.548 

6.6035937E+00 

2. 1665290E+01 

3.8216101E+01 

35.000 

114.829 

6.621H877E+00 

2.1725309E+01 

3.8031500Et-01 

36.000 

118.109 

6. 63691 83E-^00 

2. 1774622E^01 

3.7847361E401 

37.000 

121.390 

6.6492621E+00 

2 .1815120E+01 

3.7663600Ef01 

38.000 

124.671 

6. 6593956E+00 

2.1848366E+01 

3.74HO149E+01 

39.000 

127.952 

6.66771 19E+00 

2.1875651Et-01 

3.72969516+01 

40.000 

131.233 

6.6745353E+00 

2.1898037E+01 

3.7113960E+01 

41.000 

134.514 

6.6801327E+00 

2.1916401E-*'0i 

3.69311396+01 

42.000 

137.794 

6.6847235E+00 

2.1931463E+01 

3.6748456E+01 

43.000 

141.075 

6.6884R83E+U0 

2.1943815E+01 

3.6565889E+01 

44.000 

144.356 

6.6915752E*-00 

2.1953942E+01 

3.638341 36+01 

45.00(' 

147.637 

6.694l062E-^00 

2 . 19622 46E-t-01 

3.6201014E+01 

46.000 

150.918 

6.6961812E+00 

2.1969054E>01 

3.6018678E+01 

47.000 

154. 199 

6.6978822E-*00 

2 .1974634E+01 

3.5836393E+01 

48.000 

157.479 

6.6992767E+00 

2. 1979210E+01 

3.5654150E+01 

49.000 

160.760 

6.700419 7E-t-00 

2.1982960E+01 

3.5471940E+01 

50.000 

164.041 

6.7013566E-*-00 

2. 1986033E+01 

3.52807606+01 
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C. Program Listing 
I. The main program 


C LINCAK ENTRV-THAJECTORV CALCULATIONS 

C MAIN program 

1 format ISFIO.OI 


2 FORMAT I6F10.0I 


11 

FORMAT 

12 AH 

SURFACE PRESSURE 

■,F10.3 

,5H 

IMB) 1 

12 

FORMAT 

(2AH 

SURFACE density 

>,F10.3 

• 12H 

lE-OS G/CC 1 ! 

13 

FORMAT 

12 AH 

SURFACE TEMPERATURE 

-,F10.3 

,8H 

(OEG.K ) ) 

lA 

FORMAT 

I2AH 

SURFACE GRAVITY 

■ ,F10.3 

,10H 

(CM/Sa.S M 

IS 

FORMAT 

12 AH 

SPECIFIC HEAT RATIO 

■,F10.3 

) 


16 

FORMAT 

I2AH 

ADIABATIC LAPSE RATE 

■,F10.3 

,11H 

IOEG.k/KMI ! 

17 

FORMAT 

«2AH 

TROPOPAUSE altitude 

-,F10.3 

,SH 

(KM! 1 

18 

FORMAT 

I2AH 

inverse scale height 

■«F10.3 

»7h 

( 1/KMI 1 

21 

FORMAT 

(2AH0ENTRY altitude 

>,F10.3 

,5H 

IKF) ) 

22 

FORMAT 

»2AH 

ENTRY velocity 

■,F10.3 

,7H 

(KF/S!) 

23 

FORMAT 

(2 AH 

ENTRY ANGLE 

■,F10.3 

,10H 

(DEGREES! ! 

2A 

FORMAT 

1 2AH 

ballistic coefficient 

• ,F10.3 

,13H 

ISLUGS/SQ.F ! ! 

31 

FORMAT 

( 9H1 altitude, AX, 8HALTITU0E 

,BX, 




1 8HPRESSURE> nx iSHPRESSUKEt UX. SHPKESSUHEi 12X, 

2 7mOENSITV,12X,7MOENSITV) 

32 FORMAT 1 7H 1 KM I , 8X t AH( KF ) , 12X t 

1 AHIMBI « 13X ,8H|P/Sg.F| , iOXt lOHI ATMSPHRS) t9X. 

2 llHIE-05 G/CC»t5X,17M(E-05 SLUGS/CU.FII 

33 FORMAT (IHOI 

3A FORMAT (F9.3tAX,F8.3»SXt 

1 IPE14.7,5X, IPE1A.7,5X,IPE1^.7,5X, 

2 1PE14.7,5X,1PE1<>.7I 

Al FORMAT l9HlALTITU0EiAX,HHALTITUDEt8X, 

I 8MVEL0CITY,11X,8MVEL0CITV,«X,15MT|ME FROM ENTRY! 

a2 format (7H IKMI ,8X,AHIKFI . llXt 

1 6H(KM/SI t 13Xt6H(KF/S) f i 3X . 6HI SEC . ) ) 

A3 FORMAT (F9.3tAX,Fe.3i5X,lPElA.7,SX,lPElA.7,3X«lPElA.7l 
READ I»PO.OO.TO«SGiSHR,ALRtHT,SHl 
READ 2,HE .VEtEA.BC.ALT.DELTA 
PRINT 11, PO 
PRINT 12,DO 
PRINT 13, TO 
PRINT 1 a,SG 
PRINT 15,SHR 
PRINT 16,ALR 
PRINT 17, HT 
PRINT 18, SHI 
SG-SG/100. 

PRINT 21, HE 

PRINT 22, VE 

PRINT 23, EA 

PRINT 2A,BC 

HE-HE*.30A80061 

VE-VEA.30A80061 

BC>BC*157. 08682 

EA«EA*.017A5329 

C0V«2.*SG*BC*S1NF(EAI 

CDT-SH|*VE*S1NF( EAI 

PRINT 31 

PRINT 32 

PRINT 33 

H«0. 

103 CALL APDP (P0,D0,T0,SHR,ALR,HT,Sh 1,H,P,DI 
HC-H«3. 2808333 
PC«P*2.08855lA 
PA-P*. 00098692 
0C-D«1.9A03233 
PRINT 3A,H,HC,P,PC,PA,D,DC 
IF IH-ALTI 101,102,102 

101 H-H^DELTA 
GO TO 103 

102 PRIKT A1 
PRINT A2 
PRINT 33 

CALL APOP(PO,00,TO,SHR,ALR,HT,SH1 ,HE,PE,DE) 

PE»PE*100. 

AE«PE/COV 
CALL E1F(AE,EIE) 

H«0. 

203 CALL APDP(PO,00,TO,SHR,ALR,HT,SHI,H,P,0) 

P-PAIOO. 

AVT-P/tDV 
V«VE*EXPF(-AVTI 
CALL E1FIAVT,E1T) 

T»(E1T-E1E)/CDT 
HC-H«3. 2808333 
VC-VA3. 2808333 
PRINT A3,H,HC,V,VC,T 
IF Ih-ALTI 201,202,202 

201 H-H4-DELTA 
GO TO 203 

202 STOP 
END 
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2 . The eubroutinet. 


SUBAOUTINe APDPIAOtOOtTOtSHH,ALM,HT.SH| tMtPiO) 
C ATMO^PHEHIC PH(SSU«i/OCNS|TV PKOPiLE 

XPNTP*SHR/(SHA-1. I 


SUSROOTINE eiPM.tn 
C exponential integral fUNCTION 

EC«. 9772196649 


XPNTD«l./(SHR-l.| 

argt«i.*alrpht/to 

PT-PU*ARGT«*XPNTP 
PS>PT«ExPPISH|«HT) 
OTbOO«ARGT**XPNTO 
OS-DT*EXPPISH|»hTI 
IF (h-hT) I01.102fl02 

101 ARG>1.*ALR*H/T0 
PaPO*ARG««XPNTP 
0«D0*ARG»*XPNT0 

fi F TIIUIM 

102 P-PS*EXPFI-SH|*M| 
0»0S*EXPF I-SHI •Ml 

return 

ENO 


Xl«l. 

FCTRL-l. 

SUM-0. 

103 fctrl-fctrl*xi 

TERM-(X«*X| I/IFCTRL*X| I 
SUM-SUM«TERM 

IF ITEKM/SUM>1.E>0S I 101.102.102 

101 EI-EC-LOGFIXI-SUM 
RETURN 

102 Xl-Xl-1. 

GO TO 103 
ENO 
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